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We show here that many of the normal state properties of the cuprates can result from the new 
charge 2e bosonic field which we have recently (Phys. Rev. Lett. 99, 46404 (2007) and Phys. 
Rev. B 77, 014512 (2008)) shown to exist in the exact low-energy theory of a doped Mott insulator. 
In particular, the 1) mid-infrared band including the non-vanishing of the restricted f-sum rule in 
the Mott insulator, 2) the T 2 contribution to the thermal conductivity, 3) the pseudogap, 4) the 
bifurcation of the electron spectrum below the chemical potential as recently seen in angle- resolved 
photoemission, 5) insulating behaviour away from half-filling, 6) the high and low-energy kinks in 
the electron dispersion and 7) T-linear resistivity all derive from the charge 2e bosonic field. 

We also calculate the inverse dielectric function and show that it possesses a sharp quasiparticle 
peak and a broad particle-hole continuum. The sharp peak is mediated by a new charge e composite 
excitation formed from the binding of a charge 2e boson and a hole and represents a distinctly new 
prediction of this theory. It is this feature that is responsible for dynamical part of the spectral 
weight transferred across the Mott gap. We propose that electron energy loss spectroscopy at finite 
momentum and frequency can be used to probe the existence of such a sharp feature. 



Key challenges facing a theory of the normal state (at 
zero magnetic field) of the copper oxide high temperature 
superconductors include 1) T-linear resistivity^, 2) the 
pseudogap 2,3,5 , 3) absence of quasiparticles 4 , and 4) the 
mid-infrared band in the optical conductivity-^ 7 -^ ' 10 ' 11 . 
Since the parent cuprates are Mott insulators, the normal 
state properties should, in principle, be derivable from 
the corresponding low-energy theory. Proper low-energy 
theories are constructed by integrating out the degrees of 
freedom far away from the chemical potential. 

While no shortage of low-energy theories has been 
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none, until recentl y 30 ' 31 , has been based on an explicit 
integration of the degrees of freedom at high energy, 
even in the simplest model for a doped Mott insulator, 
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Here i,j label lattice sites, gij is equal to one iff i,j 
are nearest neighbours, Ci„ annihilates an electron with 
spin a on lattice site i, t is the nearest- neighbour hop- 
ping matrix element and U the energy cost when two 
electrons doubly occupy the same site. The cuprates live 
in the strongly coupled regime in which the interactions 
dominate as t ~ 0.5eV and U = 4cV. As U is the largest 
energy scale, it is appropriate to integrate over the fields 
that generate the physics on the U scale. We showe d 30 ' 31 
explicitly how to perform such an integration. The corre- 
sponding low-energy theory contains degrees of freedom, 
namely a charge 2e boson, that are 1) not in the original 
model and 2) more importantly, are not made out of the 
elemental excitations. We show here that this new de- 
gree of freedom mediates many of the anomalies in the 
normal state of the cuprates as probed by electrical and 
thermal transport as well as angular-resolved photoemis- 



sion spectroscopy (ARPES). 



I. OVERVIEW OF LOW-ENERGY THEORY 

In a series of paper o 30 ' 31 , we showed how to coarse 
grain the Hubbard model cleanly for U ^> t. We ac- 
complished this by extending the Hilbert space of the 
Hubbard model and associating with the high energy 
scale a new fermionic oscillator which is constrained. 
The coupling constant for the fermionic oscillator is U. 
We showed that once the constraint is solved, we ob- 
tained exactly the Hubbard model. However, since the 
energy scale of the new osciallator is U, the low-energy 
theory is obtained by integrating over this field. We 
showed explicitly 30 ^ that the Lagrangian in the ex- 
tended Hilbert space is quadratic in the high-energy field 
and hence all integrals can be performed exactly. Let us 
define 
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and b t = J2jhj = 0*jCj>V r tr c il _ <7 with V T = 
— V| = 1. At zero frequency, the exact low-energy Hamil- 
tonian is 



where 
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which contains a charge 2e bosonic field, <pi. 

The low-energy theory acts only in the original Hilbert 
space of the Hubbard model because we integrated out 
all operators which acted in the extended space. Con- 
sequently, it is incorrect to interpret ipi as a canonical 
boson operator with an associated Fock space. Likewise, 
we should not immediately conclude that ip^ gives rise to 
a propagating charge 2e bosonic mode, as it does not have 
canonical kinetics; at the earliest, this could be generated 
at order 0(t 3 /U 2 ) in perturbation theory. Alternatively, 
we believe that (p appears as a bound degree of freedom. 
The new composite charge e state is crucial to explain- 
ing the origin of dynamical spectral weight transfer across 
the Mott gap. As a consequence, the conserved charge is 
no longer just the electron number but rather, 
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The terms containing ipi are absent from projective 
theories and represent the fact that consistent with the 
Hubbard model, the corresponding low-energy theory is 
not diagonal in any sector containing a fixed number of 
doubly occupied sites. As we have shown elsewhere 3 ^., 
the presence of double occupancy in the low-energy the- 
ory presented above is completely compatible with the 
standard derivation of the t-J model. One simply has to 
remember that the operators appearing in the t-J model 
are transforme d 32 ' 33 ' 34 fermions which are related by a 
similarity transformation to the bare electrons appear- 
ing in the Hubbard model. If the rotated fermions are 
expressed^ 5 , in terms of the original electron operators 
in the Hubbard model, double occupancy of the bare 
fermions is reinstated. In fact, the t-J model expressed 35 
in terms of the original bare electron operators contains 
processes which mix the doubly and singly occupied sec- 
tors with a matrix element t 2 /U 

([3]). However, it is common to 

interpre t 12,13 - 14 - 15 - 16 - 17 - 18,19 - 20,21 - 22,23 - 24,25 - 26 - 27 - 28 - 29 
the t — J model as a model about physically doped holes 
in a Mott insulator. Unless the operators in the t — J 
model are transformed back to the original basis, it is 
impossible to make any connection between doped holes 
in the two theories. 

Physical processes mediated by the boson include sin- 
glet motion (third term in Eq. ([3])) and motion of double 
occupancy in the lower Hubbard band (second term in 
Eq. |(3J|). The latter has a bandwidth of t and, as we will 
show, constitutes the mid-infrared band. Since Eq. Q 
retains all the low-energy degrees of freedom, a (t/U) ex- 
pansion is warranted. To leading order in t/U, M. = Sij 
and the effective low-energy Hamiltonian simplifies to 

t 2 ■ t 2 
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FIG. 1: Hopping processes mediated by the (fii operator in 
the low-energy theory obtained by explicitly integrating out 
the high energy sector. This process is in the t-J model only 
if one tranforms the electron operators back to the orignal 
Hubbard basis, retaining the terms that change the number 
of doubly occupied sites. 
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The first two terms contain the interactions in the t — 
J model as the two-site term in b-bj is proportional to 
the unprojected spin-spin interaction. In second-order 
perturbation theory, the interaction term ip\bi mediates 
the process in Fig. (flj. In this process, hole motion over 
three sites is possible only if a doubly occupied site and a 
hole are neighbours. This process is absent if one assumes 
that the no-double occupancy condition in the t-J model 
applies to the bare electrons as well. The latter is true 
exactly at U = oo. As we will see, this process, which 
is present as a result of the bosonic degree of freedom, is 
responsible for many of the anomalous properties of the 
normal state of the cuprates. 

A. Electron Spectral Function 

We analyze the nature of the electronic excitations at 
low energy by focusing on the electron spectral function. 
Since this calculation cannot be done exactly, we estab- 
lish at the outset the basic physics that a correct calcu- 
lation should preserve. We have shown previousl y 30 ' 31 
that the electron creation operator at low energy, 

4,„ - (l-ni^c^ + V^blM^c^ 

- V^M^cj^, (6) 

contains the standard term for motion in the lower Hub- 
bard band (LHB), (1 — Ui,_o-)c| CT (rii^ a Ci a in the upper 
Hubbard band (UHB) for electron doping) with a renor- 
malization from spin fluctuations (second term) and a 
new charge e excitation, c^-^-M"- Vj ■ In the lowest order 
in t/U, our theory predicts that the new excitation cor- 
responds to a.~ a Lp\, that is, a hole bound to the charge 
2e boson. This extra charge e state mediates dynamical 
(hopping-dependent) spectral weight transfer across the 
Mott gap. 

Consequently, we predict that an electron at low en- 
ergies is in a superposition of the standard LHB state 
(modified with spin fluctuations) and a new charge e state 
which is a composite excitation. It is the presence of these 



3 



two distinct excitations that gives rise to the static (state 
counting giving rise to 2x)2S, and dynamical parts of the 
spectral weight transfer. A saddle-point analysis will se- 
lect a particular solution in which ipi is non-zero. This 
will not be consistent with the general structure of Eq. 
(JOJ) in which part of the electronic states are not fixed by 
ipi- Similarly, mean field theory in which ipi is assumed 
to condense, thereby thwarting the possibility that new 
excitations form, is also inadequate. 

The procedure we adopt is the simplest that preserves 
the potential strong interactions between the Bose and 
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where FT refers to the Fourier transform and T is the 
time-ordering operation. Further, since our focus is on 
the new physics mediated by the interaction between the 
Bose and Fermi degrees of freedom and Eq. §6§ sug- 
gests that the spin-spin interaction simply renormalizes 
the electronic states in the LHB, we will neglect the spin- 
spin term. To see what purely fermionic model underlies 
the neglect of the spin-spin term in Eq. ([5]), we inte- 
grate over (fi . The full details of how to carry out such 
an integration are detailed elsewhere^. The resultant 
Hamiltonian is not the Hubbard model but rather, 

H' = H Hubb - t -Y d b\h, (9) 
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&t — J — U model in which the spin-exchange interaction 
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where 



The effective Lagrangian can be diagonalized and written 



Fermi degrees of freedom. We treat <f>i to be spatially 
independent with no dynamics of its own. This interpre- 
tation is consistent with the fact that (pi acquires dynam- 
ics only through electron motion. Under this assumption, 
the single-particle electron Green function 

G(k,u) = -iFT{Tci(t)<S(0)), (7) 

can be calculated rigorously in the path-integral formal- 
ism as 



(8) 
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is not a free parameter but fixed to J = —t 2 /U. That 
the t — J — U model with J = —t 2 /U is equivalent to a 
tractable IR model, namely Eq. ([5]) without the spin-spin 
term, is an unexpected simplification. As Mott physics 
still pervades the t — J—U model in the vicinity of half- 
filling, our analysis should reveal the non-trivial charge 
dynamics of this model. That is, the physics we uncloak 
here is independent of the spin degrees of freedom. In 
fact, what we show is that the spin-spin interaction (as 
in the hard-projected t — J model) is at best ancillary to 
many of the normal state properties of the cuprates. 



To proceed, we will organize the calculation of G(k, oj) 
by first integrating out the fermions (holding ip fixed) 
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in terms of a set of Bogoliubov quasiparticles, 

7fe| = + cos2 Q kC* k] + sin 2 9 k C- kL 
— sin 2 fi^c^ + cos 2 9 k c- k i 



(13) 
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where cos 2 9 k = ^(1 
rj k = 2(cosfeE - cosfc^), Eq = -(2/i + £)tp*ip, E k = 



, jj-). Here, afe = 2(cos/cj; + cos/c. s ), 
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s<p*(l - §a k ) and 
when 6 = 1—n — > l — Q+2tp*tp is a renormalized 
factor which originates from the correlated hopping term 
(1 - n lS )c\ a c 0<J {l — n oS ). Starting from Eq. |12|), we 
integrate over the fermions in Eq. (|10p to obtain, 



G{k,uj) = - / [D<p*][Dcp}G(k,uj,p)exp 
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where 



G(k,u,(p) = 
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is the exact Green function corresponding to the La- 
grangian, Eq. (fT2|) . The two-pole structure of G(k, u>, ip) 
will figure prominently in the structure of the electron 
spectral function. To calculate G(k,oj), we numerically 
evaluated the remaining ip integral in Eq. (|15p . Since 
Eq. (|15j) is averaged over all values of tp, we have circum- 
vented the problem inherent in mean-field or saddle-point 
analyses. Physically, Eq. (fl"5")) serves to mix (through the 
integration over <p) all subspaces with varying number of 
double occupancies into the low-energy theory. Hence, 
it should retain the full physics inherent in the bosonic 
degree of freedom. 

The spectral function for U = lOt evaluated from Eq. 
(fT6|) and displayed in Figs. (J2J) and Fig. ([3|) exhibits four 
key features. First, regardless of the doping, there is a 
low-energy kink in the electron dispersion. The enlarged 
region in Fig. ((2^) shows, in fact, that two kinks exist. 
The low-energy kink occurs at roughly 0.2t w lOOmeV. 
By treating the mass term for the boson as a variable 
parameter, we verified that the low-energy kink is de- 
termined by the bare mass. In the effective low-energy 
theory, the bare mass is t 2 /U. This mass is indepen- 
dent of doping. Experimentally, the low-energy kink^i is 
impervious to doping. Consequently, the boson offers a 
natural explanation of this effect that is distinct from the 
phonon mechanisms which have been proposed^. 

Second, a high-energy kink appears at roughly 0.5i ss 
250meV which closely resembles the experimental kink 
at 300meV— . At sufficiently high doping (see Figs. ((3^) 
and (Jib)), the high-energy kink disappears. 

Third, experimentally, the high-energy kink is accom- 
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FIG. 2: (a) Spectral function for filling n = 0.9 along the 
nodal direction. The intensity is indicated by the color 
scheme, (b) Location of the low and high energy kinks as 
inidcated by the change in the slope of the electron disper- 
sion, (c) The energy bands that give rise to the bifurcation 
of the electron dispersion. 



panicd by a splitting of the electron dispersion into two 
branches-®. As is evident, this is precisely the behaviour 
we find below the chemical potential. The energy dif- 
ference between the two branches achieves a maximum 
at (0, 0) as is seen experimentally. A computation of 
the spectral function at U = 20t and n = 0.9 reveals 
that the dispersion as well the bifurcation still persist. 
Further, the magnitude of the splitting does not change, 
indicating that the energy scale for the bifurcation and 
the maximum energy splitting are set by t and not U. 
The origin of the two branches is captured in Fig. ([2fc). 
The two branches below the chemical potential corre- 
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FIG. 3: Spectral function for two different fillings (a) n — 0.8 
and (b) n = 0.4 along the nodal direction. The absence of 
a splitting in the electron dispersion at n — 0.4 indicates 
the bifurcation ceases beyond a critical doping. The spec- 
tral functions for two different values of the on-site repulsion, 
(c)U = lOt and (d)U = 20t for n = 0.9 reveals that the high- 
energy kink and the splitting of the electron dispersion have 
at best a weak dependence on U. This indicates that this 
physics is set by the energy scale t rather than U. 



spond to the standard band in the LHB (open square in 
Fig. on which ip vanishes and a branch on which 

ip 7^ (open circles in Fig. (J2}:). The two branches in- 
dicate that there are two local maxima in the integrand 
in Eq. (|15p . One of the maxima, ip = 0, arises from 
the extremum of G(k,ui,(p) whereas the other, the ef- 
fective free energy (exponent in Eq. (|15p ) is minimized 
(ip 7^ 0). Above the chemical potential only one branch 
survives. The split electron dispersion below the chem- 
ical potential is consistent with the composite nature of 
the electron operator dictated by Eq. ©. At low ener- 
gies, the electron is a linear superposition of two states, 
one the standard band in the LHB described by excita- 
tions of the form, c ia (l — nig) and the other a composite 
excitation consisting of a bound hole and the charge 2e 
boson, Ci&tp\. The former contributes to the static part 
of the spectral weight transfer (2x) while the new charge 
e excitation gives rise to the dynamical contribution to 
the spectral weight transfer. Because the new charge 
e state is strongly dependent on the hopping it should 
disperse as is evident from Fig. ([3]) and also confirmed 
experimentally. 

The formation of the composite excitation, Ci S p\ leads 
to a pseudogap at the chemical potential primarily be- 
cause the charge 2e boson is a local non-propagating de- 
gree of freedom. The spectral functions for n = 0.9 and 
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FIG. 4: (a)dc electrical resistivity as a function of temper- 
ature for n = 0.9 (b) Setting the bosonic degree of free- 
dom to zero kills the divergence of the resistivity as T —* 0. 
This suggests that it is the strong binding between between 
the fermionic and bosonic degrees of freedom that ultimately 
leads to the insulating behaviour in the normal state of a 
doped Mott insulator. 



n = 0.8 both show an absence of spectral weight at the 
chemical potential. Non-zero spectral weight resides at 
the chemical potential in the heavily overdoped regime, 
n = 0.4, consistent with the vanishing of the pseudo- 
gap beyond a critical doping away from half-filling. Be- 
cause the density of states vanishes at the chemical po- 
tential, we expect that the electrical resistivity to diverge 
as T — > 0. Such a divergence is shown in Fig. an d is 
consistent with our previous calculations of the dc resis- 
tivity using a local dynamical cluster method^. When 
the boson is absent (Fig. (0b)), localization ceases. Al- 
though this calculation does not constitute a proof, it 
is consistent with localisation induced by the formation 
of the bound composite excitation, Ci S p\. This state of 
affairs obtains because the boson is not an inherently dy- 
namical excitation. 

Finally, the imaginary part of the self energy at differ- 
ent temperatures is shown in Fig. (0. At low temper- 
ature (T < t 2 /U), the imaginary part of the self-energy 
at the non-interacting Fermi surface develops a peak at 
u) = 0. At T = 0, the peak leads to a divergence. This is 
consistent with the opening of a pseudogap. As we have 
pointed out earlier—, a pseudogap is properly identified 
by a zero surface (the Luttinger surface) of the single- 
particle Green function. This zero surface is expected 
to preserve the Luttinger volume if the pseudogap lacks 
particle-hole symmetry as shown in the second of the fig- 
ures in Fig. ([5]). 



B. Mid-Infrared Band 

The mid-infrared band (MIB) in the cuprates is a sur- 
prise because the optical conductivity in a doped Mott 
insulator is expected to be non-zero either at the far- 
infrared or the ultra-violet or upper-Hubbard-band scale. 
While many mechanisms have been proposed* 1 , no expla- 
nation has risen to the fore. Experimentally, the intensity 
in the MfB increases with doping at the expense of spec- 
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FIG. 5: The imaginary part of the self energy as the function 
of temperature for n = 0.7. A peak is developed at u> = 
at low temperature which is the signature of the opening of 
the pseudogap. The density of states explicitly showing the 
pseudogap is shown in adjacent figure. 



tral weight at high energy and the energy scale for the 
peak in the MIB is the hopping matrix element t. Since 
the MIB arises from the high-energy scale, the current 
theory which accurately integrates out the high energy 
degrees of freedom should capture this physics. To obtain 
a direct link between the conductivity and the spectral 
function, we work in the non-crossing approximation 



&xx{u) = 27re 2 J d 2 k J duj'(2t sin k x 



uJ 



A(ui+aj',k)A(ui',k) 
(17) 



to the Kubo formula for the conductivity where f(u>) is 
the Fermi distribution function and A(lo, k) is the spec- 
tral function. At the level of theory constructed here, 
the vertex corrections are all due to the interactions with 
the bosonic degrees of freedom. Since the boson acquires 
dynamics only through electron motion and the leading 
such term is 0(t 3 /U 2 ), the treatment here should suffice 
to provide the leading behaviour of the optical conduc- 
tivity. 

The optical conductivity shown in Fig. ([6]) peaks at 
w/fsi JSt forming the MIB. As the inset indicates, w max 
is an increasing function of electron filling (n) whereas 
the integrated weight 
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is a decreasing function. However, N e g does not vanish 
at half-filling indicating that the mechanism that causes 
the mid-IR is evident even in the Mott state. Here we set 
the integration cutoff to f2 c = 2t = 1/ m* . Both the mag- 
nitude of f2 max and its doping dependence as well as the 
electron filling dependence of the integrated weight are 
consistent with that of the mid-infrared band in the op- 
tical conductivity in the cuprates ^ 9 ' 10 ! 11 . To determine 
what sets the scale for the MIB, we studied its evolution 




FIG. 6: Optical conductivity as a function of electron fill- 
ing, n. The peak in the optical conductivity represents the 
mid-infrared band. Its origin is mobile double occupancy in 
the lower-Hubbard band. The insets show that the energy 
at which the MIB acquires its maximum value, oj max is an 
increasing function of electron filling. Conversely, the inte- 
grated weight of the MIB decreases as the filling increases. 
This decrease is compensated with an increased weight at 
high (upper-Hubbard band) energy scale. 



as a function of U. Figure ([7]) verifies that w max is set 
essentially by the hopping matrix element t and depends 
only weakly on J. The physical processes that determine 
this physics are determined by the coupled boson-Fermi 
terms in the low-energy theory. The ipjc^Cii term has a 
coupling constant of t whereas the ip\bi scales as t 2 /U. 
Together, both terms give rise to a MIB band that scales 
as w max /i = 0.8 - 2.2U/U (see inset of Fig. ©). Since 
t/U w O(.l) for the cuprates, the first term dominates 
and the MIB is determined predominantly by the hop- 
ping matrix element t. Within the interpretation that 
tp represents a bound state between a doubly occupied 
site and a hole, second order perturbation theory with 
the (f\bi term mediates the process shown in Fig. (TT]). It 
is the resonance between these two states that results in 
the mid-IR band. Interestingly, this resonance persists 
even at half-filling and hence the non-vanishing of N c s 
at half-filling is not evidence that the cuprates are not 
doped Mott insulators as has been recently claimed 4 ^. 
In their work, Comanac and colleagues 4 ^ used a single- 
site dynamical mean-field approach. In such approaches, 
near-neighbour correlations are absent. 

Since the physics in Fig. (TT]) is not present in projective 
models which prohibit double occupancy in the Hubbard 
basis (not simply the transformed fermion basis of the 
t-J model), it is instructive to see what calculations of 
the optical conductivity in the t — J model reveal. All 
existing calculation s 10 ! 42 ' 43 on the t — J model find that 
the MIB scales as J. In some of these calculations, su- 
perconductivity is needed to induce an MIB 43 also at an 
energy scale of J. Experimentally&ili 1 -, it is clear that 
the MIB is set by the t scale rather than J. In fact, 
since the MIB grows at the expense of spectral weight 
in the upper-Hubbard band, it is not surprising that the 
t-J model cannot describe this physics as first pointed 



FIG. 7: Evolution of the optical conductivity for n = 0.9 as 
U is varied. The inset shows the functional form that best 
describes cj max . The dominant energy scale is the hopping 
matrix element t since t/U for the cuprates is 0(1/10). 

out by Uchida, et ali The physical mechanism we have 
identified here, Fig. |T]) clearly derives from the high en- 
ergy scale, has the correct energy depenence, and hence 
satisfies the key experimental constraints on the origin of 
the MIB. Since the physics in Fig. dXJ) is crucial to the 
mid-IR, it is not surprising that single-site analyses^ fail 
to obtain a non-zero intercept in the extreme Mott limit. 
The non-zero intercept of N e g is a consequence of Mot- 
tness and appears to be seen experimentally in a wide 
range of cuprates£ > 44 ' 45 i 46 

C. Dielectric function: Experimental Prediction 

In the previous section, we have calculated the elec- 
tronic spectral function which shows that there are two 
branches below the chemical potential. Such physics is 
explained by the formation of a new composite excita- 
tion, representing a bound state, consisting of a bound 
hole and a charge 2e boson, tp\ci S . We demonstrated that 
for the MIB in the optical conductivity such an excita- 
tion also appears. In principle, these composite charge 
excitations should show up in all electric response func- 
tions, for example, the energy loss function, 31/e(w,q), 
where e(ui, q) is the dielectric function. We show here 
that this is the case. 

To this end, we calculate the inverse dielectric function, 

e(w,q) v p J 

x A(u+u',p + q)A(u;',p), (19) 

using the non-crossing approximation discussed earlier. 
Our results are shown in Fig.® for n = 0.9 and n = 0.6 
for q along the diagonal. Two features are distinct. First, 
there is a broad band (red arrow in Fig. (|5J|) with the 
width of order t that disperses with q for both doping 




FIG. 8: The dielectric function, — Syl/e(o», q) for q along the 
diagonal direction is shown for (a) n = 0.9 and (b) n = 0.6. 
Note only the broad feature indicated by the red arrow at 
n — 0.9 persists at n = 0.6. 



levels. It is simply the particle-hole continuum which 
arises from the renormalized bare electron band. The 
band width is doping dependent as a result of the renor- 
malization of the band with doping. More strikingly, for 
n = 0.9, a sharp peak exists at w/f w .It. It disperses 
with q, terminating when q — ► (w, w). Physically, the 
sharp peak represents a quasiparticle excitation of the 
composite object, (p\cia, the charge 2e boson and a hole. 
Therefore, we predict that if this new composite charge 
excitation, <p\ci S , is a real physical entity, as it seems 
to be, it will give rise to a sharp peak in addition to the 
particle-hole continuum in the inverse dielectric function. 
Since this function has not been measured at present, our 
work here represents a prediction. Electron-energy loss 
spectroscopy can be used to measure the inverse dielec- 
tric function. Our key prediction is that momentum- 
dependent scattering should reveal a sharp peak that ap- 
pears at low energy in a doped Mott insulator. We have 
checked numerically the weight under the peaks in the 
inverse dielectric function and the sharp peak is impor- 
tant. Hence, the new charge e particle we have identified 
here should be experimentally observable. 

The two dispersing particle-hole features found here 
are distinct from a similar feature in stripe models^. In 
such models the second branch^ has vanishing weight 
and whereas in the current theory both features are of 
unit weight. 
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FIG. 9: (a) Heat capacity, CV, and (b) thermal conductivity, 
t, calculated at n = 0.9. The solid lines are a fit to T 2 . Insert: 
Density of states for U = Wt axe evaluated at n = 0.9 and 
n = 0.7 respectively. 



D. Heat conductivity and heat capacity 

As shown by Loram and collaborators 4 ^, the heat ca- 
pacity in the cuprates in the normal state scales as T 2 . 
Quite generally, a density of states that vanishes linearly 
with energy, that is, a V-shaped gap, yields a heat ca- 
pacity that scales quadratically with temperature. The 
coefficient of the T 2 term is determined by the slope of 
the density of states in the vicinity of the chemical poten- 
tial. The magnitude of the T 2 term should diminish as 
the doping increases because the slope of the density of 
states decreases as the pseudogap closes. As we showed 
in the previous section, the boson creates a pseudogap. 
The energy dependence of the gap is shown in the inset 
of Fig. ©. Quite evident is the linear dependence on 
energy. The resultant heat capacity shown in Fig.([9^), 
calculated via the relationship C v = , where the inter- 



nal energy, E, is 



E 



du)D(uj)uj f{u>) 



(20) 



and D(uj) = ^ k A(u>, k) is the density of states displays a 
perfectly quadratic temperature dependence in the dop- 
ing regime where the pseudogap is present as is seen 
experimentally 48 .. As it is the boson that underlies the 
pseudogap, it is the efficient cause of the T 2 dependence 
of the heat capacity. In our theory, the steeper slope 
occurs at smaller doping which gives rise to the largest 
heat capacity at half filling. This doping dependence of 
the heat capacity seems to contradict the experimental 
observations^,. A key in determining the magnitude of 
the heat capacity is the spin degrees of freedom. As we 
have focused entirely on the bosonic degree of freedom 
and not on the contribution from the spin-spin interac- 
tion terms, we have over-estimated the kinetic energy. 
Such terms, though they do not affect the pseudogap 
found here (from Eq. © it is clear that the spin-spin 
terms renormalize the standard fermionic branch in the 
lower-Hubbard band leaving the new state mediated by 
ifii untouched), do alter the doping dependence 4 ^. 

Additionally, the thermal conductivity, t(T), can be 



calculated using the Kubo formula in non-crossing ap- 
proximation, 



r(T) = 



4fc B T 



The thermal conductivity shown in Fig.@ scales as T 2 
which is identical to that of the heat capacity. However, 
the system exhibits a larger thermal conductivity as the 
doping increases in contrast to the heat capacity which 
is decreasing as the doping increases. Physically, this 
signifies that the carriers are more mobile as the doping 
increases. 



E. T-linear Resistivity 

Over a funnel-shaped region in the T — x plane, the 
resistivity displays the anomalous linear-T dependence. 
The standard explanation 4 ^ attributes T— linear resisi- 
tivity to quantum criticality. However, one of us has re- 
cently shown 4 ^ that under three general assumptions, 1) 
one-paramater scaling, 2) the critical degrees of freedom 
carry the current and 3) charge is conserved, the resis- 
tivity in the quantum critical regime takes the universal 
form, 



O 2 

a(uj = 0) = 4" 
n 



k B T 

he 



(d-2)/z 



(21) 



Consequently, T— linear resistivity obtains (for d=3) only 
if the dynamical exponent satisfies the unphysical con- 
straint z < 0. The inability of Eq. (|2"Tj) to lead to a 
consistent account of T— linear resistivity signifies that 
1) either T— linear resistivity is not due to quantum crit- 
icality, 2) additional non-critical degrees of freedom are 
necessarily the charge carriers, or 3) perhaps some new 
theory of quantum criticality can be constructed in which 
the single-correlation length hypothesis is relaxed. 

We show that the low-energy theory presented here 
contains elements of both (2) and (3) which lead to 
T— linear resistivity. The formation of the pseudogap 
and the divergence of the electrical resistivity are highly 
suggestive that a bound state between a hole and the 
charge 2e bosonic field, namely the tp^Cig particle dis- 
cussed earlier in the electron operator and the new fea- 
ture in the dielectric function. Let assume this state of 
affairs obtains and the binding energy is Eb ■ As a bound 
state, Eb < 0, where energies are measured relative to 
the chemical potential. Upon increased hole doping, the 
chemical potential decreases. Beyond a critical doping, 
the chemical potential, crosses the energy of the bound 
state. At the critical value of the doping where Eb = 0, 
the energy to excite a boson vanishes. The critical region 
is dominated by electron-boson scattering. In metals, it 
is well-known^ that above the Debye temperature, the 
resistivity arising from electron-phonon scattering is lin- 
ear in temperature. In the critical regime, the current 
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FIG. 10: Proposed phase diagram for the binding of the holes 
and bosons that result in the formation of the pseudogap 
phase. Once the binding energy vanishes, the energy to ex- 
cite a boson vanishes. In the critical regime, the dominant 
scattering mechanism is still due to the interaction with the 
boson. T- linear resistivity results anytime T > u>b, where ujb 
is the energy to excite a boson. To the right of the quantum 
critical regime (QCP), the boson is irrelevant and scattering 
is dominated by electron-electron interactions indicative of a 
Fermi liquid. The QCP signifies the end of the binding of 
fermi and bosonic degrees of freedom that result in the pseu- 
dogap phase. 



mechanism yields a T— linear resisitivity for the same rea- 
son. Namely, in the critical region, the energy to create 
a boson vanishes as shown in Fig. (fTUj) and hence the 
resistivity arising from electron-boson scattering should 
be linear in temperature. This mechanism is robust as it 
relies solely on the vanishing of the boson energy at criti- 
cality and not on the form of the coupling. To the right of 
the quantum critical point, standard electron-electron in- 
teractions dominate and Fermi liquid behaviour obtains. 
In this scenario, the quantum critical point coincides with 
the termination of the pseudogap phase, or equivalently 
with the unbinding of the bosonic degrees of freedom. 



F. Final Remarks 

We have shown here that 1) T-linear resistivity, 2) high 
and low-energy kinks in the electron dispersion, 3) pseu- 



dogap phenomena, 4) mid-infrared absorption, 5) T 2 heat 
capacity, and 6) the bifurcation of the electron addition 
spectrum all emerge from the bosonic degree of freedom 
which exists in the exact low-energy theory of a doped 
Mott insulator. All the features found here arise from 
the charge rather than from the spin dynamics and hence 
have nothing to do with the spin-spin interaction of the 
t — J model. The essential feature of the charge 2e boson 
is that it gives the electron substructure, as is observed 
experimentally. Further, it does so entirely from a collec- 
tive mode arising from the strong correlations. No lattice 
phenomena need be invoked. Consequently, we attribute 
the anomalies of the normal state of the cuprates to a 
purely electronic cause which is captured by the physics 
of the exact-low energy theory of a doped Mott insula- 
tor. The key experimental prediction which is sufficient 
to falsify this theory is the existence of a sharp quasiparti- 
cle excitation in the inverse dielectric function. The new 
quasiparticle represents the collective motion of a dou- 
ble occupancy and a hole as depicted in Fig. |T]) which 
emerges naturally in the exact low-energy theory as a 
charge 2e boson bound to a hole. The existence of such a 
quasiparticle underscores the fact that double occupancy 
at low energy requires a hole in its immediate vicinity. 
Since this prediction is sharp, experimental falsification 
of this theory should be straightforward. 
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